Sum rule for a Schiff-like dipole moment 
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The energy-weighted sum rule for an electric dipole transition operator of a Schiff type differs 
from the Thomas-Reiche-Kuhn sum rule by several corrective terms which depend on the number 
of system components, A/". The deviations are evaluated by two distinct approaches. One of them is 
based on the charge density expansion around an uniform distribution and provides an analytical A/" 
dependence for the sum rule while, the other one hinges of the RPA approach and yields compact 
expressions for corrections as function of the RPA amplitudes. Although the formalism might be 
considered for any many body system, we applied it for illustration, to the case of Na clusters. One 
concludes that the RPA results for Na clusters obey the modified TRK sum rule. 
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The Thomas-Reiche-Kuhn (TRK) sum rule [H, [3, Q 
has been widely used in various contexts of electron exci- 
tations in atoms, molecules, and solids. Indeed^he rela- 
tion to the photoabsorbtion oscillator strength [1, [E S 13 
makes it quite useful in interpreting various collective fea- 
tures of the dipole states. 

The many body properties like collective excitations in 
atoms, molecules or clusters are studied using different 
collective models. Thus, the random phase approxima- 
tion (RPA) formalism has been extensively used and ex- 
cellent results for the photoabsorbtion cross section spec- 
tra, especially in the atomic clusters domain [1, [^, [loj 
were obtained. 

In a previous publication [l^ , some of the many body 
features of the small and medium sodium clusters were 
studied within the RPA approach using the projected 
spherical single-particle basis defined in Ref.[ll|. The 
RPA wave functions were used to treat the dipole tran- 
sitions which led to the photoabsorbtion cross section 
spectra. Also, the system static electric polarizability 
was analytically determined. The nice virtue of the RPA 
approach with spherical single particle basis consists of 
that it satisfies TRK sum rule. 

This is however true for the electric dipole moment, 
which is not the case in the above quoted paper where, 
indeed, a modified dipole operator, similar with Schiff- 
like dipole moment [ij, [H, El, [Si is used. The correc- 
tions to the dipole operator accounts for the screening 
effect, caused by the electronic shells, on the motion of 
the valence electrons in the mean field determined by the 
nuclear charge. Moreover the approach of Ref.^] uses 
a projected spherical single particle basis which allows 
for an unified description of spherical and deformed clus- 
ters. Here we address the question whether these specific 
features require a modification of the TRK sum rule. 

Thus, in the present paper we derive a new sum rule 
which corresponds to the Schiff-like dipole moment. 

Within the RPA formalism, for any Hermitian operator 
M, the following sum rule holds 



where i?„ are the RPA energies associated to the many 
body Hamiltonian H. Here the state |0) is the RPA 
phonon vacuum, while |1„) denotes the single phonon 
state |1„) — Cl\0). We are interested in those features 
of atomic clusters which are determined by the valence 
electrons. These move in a mean field, determined by 
the ionic core, and interact among themselves through 
a Coulomb force. The residual two-body interaction 
(16, 17] can be expanded in multipole series from which 
only the dipole term is relevant and therefore consid- 
ered. The many body Hamiltonian describing the system 
is treated within the RPA formalism which defines 
the phonon operator: 



ct(i,^) = J2 [x^^clc^h^ - r;uclcp)i^' 
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by means of the equations: 

[Cn, Cj'i'] — Snn' ■ 



(2) 
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Here 



Xpf^ and Y^f^ are the RPA forward-going and 



backward-going amplitudes corresponding to the n-th or- 
der solution of the RPA equations. The reduced proba- 
bility for the dipole transition |0) |1+) can be written 
in terms of the RPA phonon amplitudes and the ph ma- 
trix elements of the transition operator : 



where 



s(i?i,o+^i+) = |(o||A^(ii;i)||i+)| 



(0||A^(i?l)||l+) = J2i,{p\mEl)\\h) 



(5) 
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are the reduced matrix elements of the dipole operator 
M{E1), between the specified RPA states. 



^(£;„-^o) (0|Af|l„ 



' = i(0|[[M,i7],M]|0), (1) 



[1] Throughout this paper the Rose's convention for the reduced 
matrix elements are used. 
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Instead of the usual transition dipqle operator, a SchifF- 
like moment operator [l^, [H, [lB| was used. 



(7) 



Here rg is the Wigner-Seitz radius and have the value 
of 3.93 a.u. for Na clusters. The corrective component, 
involved in the dipole operator, relates particle and hole 
states characterized by AiV — 3, which results in mod- 
ifying the strength distribution among the RPA states. 
Such an effect would be however obtained even for the 
dipole transition operator, if the mean field potential for 
the single particle motion involves higher powers of the 
radial coordinate. 

Reckoning the double commutator from Eq. ([T]), cor- 
responding to the transition operator ([7]) one obtains: 
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The terms correcting the standard TRK sum rule are 
the expected values of the radius powers and , 
in the RPA ground state. Obviously, these corrective 
terms induce an additional M dependence for the energy- 
weighted sum rule. As for the left hand side of Eq.®, 
this can be directly calculated using the RPA output 
data, like energies and transition probabilities The 
terms from the right hand side of Eq.® are evaluated 
alternatively through two distinct methods, which are 
briefly presented below. 

A ) The boson expansion method. The correction terms 
from the right hand side of the equation ([S]) , which will be 
hereafter denoted as S{Af), can be expressed in terms of 
particle-particle (pp) and hole-hole (hh) transition com- 
ponents. Indeed, the ph transition components give van- 
ishing contributions when they are averaged with the 
RPA ground state. Therefore in the second quantiza- 
tion, the needed one-body operator f™ can be written 
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{pn)cicp+j2{h\nh)cicH. (9) 
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Following the boson expansion procedure, the fermion 
density operators c^Cp and c\ch can be expressed in terms 
of the RPA phonon operators Cj^(C„). Thus, the particle 
and hole density operators can be written as 

= E <CnCl clcn = E b':,C^Cl (10) 

n n 

where the coefficients and h^^ have the expressions: 

< = (0|[[Ct,ctcJ,C„]|0), h^^ = (0|[[Ct,ctcJ,C„]|0). 

(11) 



In this way the modified dipole sum rule becomes: 
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where the following notation has been used: 

r\ph)^{p\r''\p)-{h\r'^\h),k^2A. (13) 

As we have already mentioned, a specific ingredient of 
the approach from Ref. [l3| is the use of the projected 
spherical single particle basis in the RPA formalism. 

B) Electron density approach. Static electric polariz- 
abilities results based on the calculations for the number 
of spilled out electrons [l3|, agree quite well with the 
corresponding experimental data. These spilled out elec- 
trons produce a screening effect against external fields 
which results in changing the classical result for the po- 
larizability. The basic assumption in accounting the 
spilled out electrons is the fact that the electron den- 
sity is not going sharply to zero at the cluster surface, 
but is gradually decreasing and moreover extends signifi- 
cantly beyond the jellium edge. The same argument can 
be brought for the correction terms of the S{N') con- 
taining averages of the radius powers with the RPA vac- 
uum states. The electronic density has a constant cen- 
tral part, enfolded by a diffuse region, of a width equal 
to a, where the electron density tends smoothly to zero. 
Guided on some parallelism between atomic clusters and 
nuclear systems [H, [l^ , the average of r™ with the RPA 
vacuum state is approximated by folding r™ with a lo- 
calization probability density with spherical symmetry, of 
Fermi distribution type, having a profile similar to that 
of the charge density, mentioned above. In this way the 
average may be written in the form of a power series in 
the variable f [lO]: 



m + 3 
771 = 2,4. 
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Here R — rgM^^^ is the radius of the cluster with TV 
atoms, with being the Wigner-Seitz radius, which for 
the Na clusters has the value of 3.93 a.u.. a is a param- 
eter defining the thickness of the diffusion region. 

In this way the sum S{N) can be written in the fol- 
lowing way 



S{N) 
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One notices that besides the number of atoms depen- 
dency, this expression involves only universal constants, 
which makes Eq. (jl5p to be, indeed, a real sum rule. 
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Contrary to the case of nuclear systems, where the 
thickness of the diffusion region is approximately the 
same for all nuclei, for atomic clusters the parameter a is 
expected to have a TV dependency due to the long range 
character of the two body interaction. This dependence 
was determined by interpolating the values of a satisfying 
Eq. (HID for 8 < A/" < 40. The resuh is: 

a{N) = -0.975157 - 0.01121387V^/^ + 0.3605187V^/^ 

(16) 

Inserting the expression pB]) of a(N) in Eg. p^ . one ob- 
tains the final expression for the sum rule. 

The sum rule is always a serious test for any model 
calculation and in particular for the RPA formalism of 
Ref.[l^ which uses a projected spherical single particle 
basis, appropriate for an unified description of spheri- 
cal and deformed clusters. The energy weighted sum of 
the dipole transition probabilities was calculated with the 
RPA energies and matrix elements for the transition op- 
erator. The right hand side of the sum rule equation, 
called S{M) was alternatively calculated with the expres- 
sions (12) and (15), respectively. The two sets of results 
corresponding to the mentioned options for S(M), are 
plotted in Figs. 1 and 2, respectively. 

The RPA calculation of S{M) does not yield an explicit 
TV dependence. However, this dependence is involved in 
the many body formalism by means of the Fermi en- 
ergy, the oscillator length, the matrix elements in the 
projected single particle basis, the space of the single 
particle states used in the RPA calculation. Note that 
Fig. 1 shows a good agreement for medium clusters. The 
deviation is significant for large as well as for small clus- 
ters although both curves exhibit similar pattern con- 
cerning the oscillating behavior. The discrepancies may 
indicate that higher order boson expansion terms are nec- 
essary in order to improve the agreement. However, these 
terms would bring a certain inconsistency to the formal- 
ism since the energy weighted sum is evaluated within the 
RPA approach and therefore is a quadratic expression of 
the dipole phonon amplitudes, while higher order terms 
of S{J\f) would be higher order polynomials of the men- 
tioned amplitudes. Moreover, in order that Eq.(12) plays 
the role of a sum rule it is necessary that S{N') exhibits a 
model independent expression. On the other hand adding 
higher order terms to 5 (TV) increases the specific fea- 
ture of the chosen RPA approach. In contradistinction 
to these aspects, the second procedure described above, 
yields indeed, a model independent expression for S'(TV). 
The results for this situation are shown in Fig. 2. 

Fig. 2 indicates a very good agreement between 5(TV) 
and the energy weighted sum. Actually this is the main 
result of our investigation. The sum rule for the modified 
dipole operator is equal to 5(TV) given by Eq.(15) with 
the diffusion parameter a from (16). One notices that 
S'(TV) is a six order polynomial in TV^^'^. The second 
remark refers to the fact that the new sum rule is different 
from the TRK sum rule, which is linear in TV. 

Now we address the question whether the experimental 
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FIG. 1: The calculated energy- weighted sums of the reduced 
dipole transition probabilities (open circles) for Na clusters 
with 8-40 atoms are compared with the numerical results of 
S(TV) given by the RPA approach (black triangles). 
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FIG. 2: The calculated energy- weighted sums of the reduced 
dipole transition probabilities (black circles) for Na clusters 
with 8-40 atoms are compared with the numerical results of 
S'(TV) given by Eq. (15) and represented by a solid line. 



data for the energy weighed sum (EWS) agree with the 
calculated values. Of course the theoretical dipole EWS 
depends on the dimension of the single particle space in- 
volved in the calculation. If the single particle states have 
a complex structure, as for example happens in Ref.jl], 
then the states with ATV = 3 are connected by the stan- 
dard dipole operator and consequently a large number of 
major shells are necessary in order to saturate the sum 
rule. If a projected spherical single particle basis is used. 
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ical weighted sum calculated within the RPA approach, 
which uses a projected spherical single particle basis p^ . 



FIG. 3: The calculated (Th.) energy-weighted sum (EWS) of 
the SchifF-like dipole transition probabilities (open circle) for 
A'^o clusters are compared with the corresponding experimen- 
tal data (cross) extracted from the available data concerning 
the photoabsorbtion cross section. 



as we did in Ref . [To| . the dipole matrix elements with 
AA^ = 3 are vanishing. In order to account for the dipole 
strength brought by the non-vanishing AN = 3 matrix 
elements we use an effective dipole operator which is the 
Schiff-like dipole operator. Therefore, the Schiff dipole 
moment is what we propose to describe the experimental 
dipole EWS. In order to complete the picture we have to 
know the experimental value for the EWS rule and check 
whether this confirms that the found sum rule is valid. 
Actually, this can be extracted form the experimental 
photoabsorbtion cross section. Indeed, interpolating the 
discrete values of the experimental photoabsorbtion cross 
sections, given as function of the excitation energy, by 
a smooth curve and integrating the result with respect 
to the energy in the interval [0, oo), one obtains the area 
A{M) ■ The experimental EWS is proportional to A{Af) . 
We assume that the JV depending proportionality factor 
is the same as for the calculated EWS, i.e. T{M). This 
assumption is grounded by the fact that T{N) is a model 
independent quantity and moreover assures similar nor- 
malization for the total cross section as in the schematic 
calculation. Thus the experimental EWS is defined by: 



[EWS]^,^^^^gn^)A{M). 



(17) 



The quantity A{Af), extracted from the data of Refs.|26l. 
[131, varies between 0.582 (eV)!^ for A/" = 9 and 0.387 
(eV)A^, for M = 19. Here the constant factor is ^ = 
70.1591912[e^]. This value was obtained by equating the 
two expressions of ESW , for one chosen cluster. Thus, 
the constant Q yields a normalization of EWS which ac- 
counts for the " missing sum rule" , noticed experimentally 
[2^. In Fig. 3 this quantity is compared with the theoret- 



[EWS]^^^ = 



EnB{Q^ 



(18) 
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The agreement between the two EWS^s, shown in Fig. 3, 
is a guarantee that [EWS]^.^^^ obeys, indeed, the equa- 
tion: 



[EWS]^^^=S{N). 



(19) 



Finally, we want to mention that the octupole correction 
to the dipole transition operator was also used in con- 
nection with the description of the electric dipole transi- 
tions in nuclear systems. Thus, in Ref.[2l| it is pointed 
out that adding the octupole correction, the agreement 
with experimental data concer ning the El transitions is 
substantially improved. In Ref. [22] a similar effect is ob- 
tained by modifying the many body wave function due to 
the octupole interaction which is considered in addition 
to an isovector-dipole interaction. The new components 
are connected by the standard dipole operator and con- 
sequently modifies the El transition rates. Of course, 
an energy-weighted sum rule associated to the dipole 
transition operator holds also for nuclear systems. Ob- 
viously,_changing the transition operator, as it happened 
in Ref. 21], the corresponding sum rule which should be 
valid is the one obtained in the present letter. 

Another interesting example is the N — Z sum rule, 
which holds for the single beta transition. This rule 
says that for a single beta decaying nucleus, the differ- 
ence between the /3~ and /3+ strengths should be equal 
to 3(A^ — Z). This sum rule is true, for example, for 
Gamow- Teller (GT) dipole transitions and is exactly sat- 
isfied within the proton-neutron quasiparticle RPA ap- 
proach. Extensions of the microscopic formalisms to the 
double beta decay 2vf3(3, showed that in order to describe 
the transition rates, it is necessary to improve the wave 
functions of the mother nucleus as well as the GT dipole 
states, by adding anharmonic effects. However, these cor- 
rections violate drastically the N — Z sum rule. In the 
spirit of the present paper, we open the question whether 
the Gamow- Teller proton-neutron interaction could be 
extended by adding an octupole component such that to 
the new proton-neutron interaction a modified N ^ Z 
sum rule corresponds. This would make the inclusion of 
anharmonic effects which, as a matter of fact, violates 
the Pauli principle, unnecessary. 

The final conclusion is that the Schiff-like dipole mo- 
ment used for the RPA description of the photoabsorb- 
tion cross section spectrum, satisfies an extended TRK 
dipole sum rule. The saturation of the extended sum rule 
is a positive test for the single particle basis as well as for 
the dimension of the dipole ph space involved in the RPA 
description. Also we showed that the experimental value 
of EWS is reproduced to a high accuracy by the sum rule 
S{JV). We pointed out that the extended TRK sum rule 
is of a general interest, being applicable also for other 
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many body systems correlated by a SchifF-like two body 
interaction. The fact that sum rules are of paramount 
importance in exploring the many body properties mir- 
rored by the multipole electric, or magnetic, transitions 
have been stressed by many authors [2^ [2J, [2^ . Here 
we showed that a sum rule may hold also for a multipole 
mixed transition operator. 
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